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We  consider  optimization  of  schedules  for  maintenance  or  repairs, 
in  order  to  minimize  long  term  average  operating  cost  or  to  maximize 
availability.  The  novelty  here  is  that  the  failure  rate  after  a  main¬ 
tenance  operation  is  a  function  of  the  system's  previously  expended 
lifetime.  This  generalizes  earlier  work  by  others  on  the  simpler  case 
where  the  future  rate  depends  only  on  the  number  of  previous  re¬ 
pairs,  but  not  on  the  times  when  they  took  place.  The  underlying 
lifetime  distributions  are  assumed  to  have  the  Weibull  form  and  two 
classes  of  maintenance  strategies  are  considered. 

The  first  case  optimizes  a  set  of  successive  maintenance  intervals 

Tt,  T2 . T„,  and  the  number  N,  where  a  replacement  by  a  new 

w  N 

system  is  made •  at  tN  =  I  T..  For  the  specified  model,  we  show 
;  1  i=l  1 

that  the  optimal  times  T.  exists  and  are  ordered  as  Tx  >  T2  >  ...  >  TN. 
Additional  properties  of  the  optimal  solutions  are  proved,  and  others 
conjectured . 

In  the  second  case,  the  period  of  periodic  maintenance  is  opti¬ 
mized  numerically.  The  main  contribution  in  that  case  is  the  formula¬ 
tion  of  the  new  failure  rate  model,  and  the  efficient  organization  of 
the  optimization  calculations. 

*  This  work  was  partially  supported  by  Grant  No.  N00014-75-C-0858 
from  the  Office  of  Naval  Research. 
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1.  INTRODUCTION 


This  work  extends  the  maintenance  models  and  policy  optimiza¬ 
tions  described  in  the  references.  In  particular,  the  ideas  in  Nguyen 
and  Murthy  (to  appear)  and  Goldstein  (1980)  were  strong  influences. 

We  consider  optimization  of  the  schedule  of  times  at  which  pre¬ 
ventive  maintenance  will  be  carried  out  on  a  system  subject  to  stochas¬ 
tic  failures.  Several  different  criteria  are  used  to  judge  the  quality 
of  a  schedule,  including  long-term  average  operating  costs  per  unit 
time,  and  various  definitions  of  availability.  The  effect  of  mainte¬ 
nance  is  represented  as  a  reduction  in  the  subsequent  failure  rate. 
One  of  our  contributions  is  an  approach  for  incorporating  the  influ¬ 
ence  of  past  operating  time  on  the  new  railure  rate  functions  which 
applies  after  a  maintenance  operation. 

A  component  or  system  with  a  failure  time  t  (random  variable) 
can  be  described  by  its  failure  rate  or  hazard  rate  r(t): 

r(t)dt  =  P[t  <  t  <  t  +  dt  |  t  >  t]  (1) 

where  the  cumulative  distribution  function  is  related  to  r(t)  as 

t 

F(t)  =  1  -  exp(-  /  r(z)  dz]  (2) 

0 


Figure  1  shows  two  forms  used  later  for  the  reduction  in  failure  rate 


at  maintenance  times  t..  In  case  a),  maintenance  reduces  r(t.)  to 


zero,  but  subsequent  growth  rate  is  greater  than  before  the  mainte¬ 
nance.  That  growth  rate  will  be  a  function  of  the  previous  operating 
time  t^.  In  case  b),  the  growth  rate  of  r(t)  is  constant,  but  mainte^, 
nance  reduces  r  at  t^  by  an  amount  dependent  on  tj. 


T 


1 
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We  define  r^(t)  between  maintenance  times  as: 

r(t)  =  r.(t  -  t._1)  t.  <  t  <  t.  .  (3) 

We  consider  r.(t)  which  are  increasing  in  t,  with  numerical  results 
based  on  the  Weibull  forms 

r^t)  =  \.ta  (4) 

or 

r.(t)  =  \t(t  +  t.^)01  -  A. 

(equivalent  to  r(t)  =  \ita  -  Aj) 

These  r.(t)  are  also  increasing  functions  of  the  maintenance  intervals 
Tj  =  (tj  "  tj-l)  for  j  <  i* 

In  some  models  we  permit  minor  repairs  when  the  system  fails. 
After  such  a  repair,  the  failure  rate  r(t)  evolves  as  if  no  failure  had 
occurred.  This  model  is  often  used  for  complex  equipment  in  which 
the  next  failure  is  unlikely  to  be  related  to  the  component  which  was 
just  repaired.  From  this  point  of  view,  the  number  of  failures  and 
minor  repairs  between  maintenance  times,  tp,  has 

'i 

EWtj.p  tp]  =  /  r.(t)  dt  (5) 

Vi 

For  both  cases  in  Fig.  1,  the  system  wears  out  (r(t)  -»  ®  as 
t  -*  »)  so  complete  replacement  (renewal)  is  eventually  required. 

Section  2  considers  a  problem  where  minor  repairs  are  permitted 
between  maintenance  times,  and  costs  are  assigned  to  renewal,  minor 
repair  and  maintenance  operations.  With  replacement  at  t^,  we  seek- 
the  best  N  and  tx,  t2,...,tN_1  to  minimize  the  average  cost  per  unit 
time  when  renewal  cycles  are  repeated  indefinitely.  This  is  mathe- 
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maticaUy  equivalent  to  maximizing  the  steady-state  availability,  by 
suitable  reinterpretation  of  the  costs. 

Section  3  assumes  that  each  failure  requires  a  complete  renewal, 
and  that  the  maintenance  times  are  equally  spaced  (often  the  case  due 
to  work  schedules).  Cycle  availability  Ac  is  defined  as  the  fraction  of 
time  that  the  system  is  up  (not  in  maintenance  or  replacement  opera¬ 
tions)  during  one  renewal  cycle.  The  maintenance  period  is  sought  to 
maximize 

P[A  >  6]  ,  1  <  5  <  0  (6) 

v 

2.  SCHEDULES  WITH  MINOR  REPAIRS 

We  assume  that  renewal  occurs  at  tN  with  a  cost  CR;  mainte¬ 
nance  occurs  at  Tx,  T2,...,  TN-1  at  a  cost  of  CM  each  time;  and 
minor  repairs  occur  at  each  failure  with  a  cost  Cr.  The  expected 
cost  per  renewal  cycle  is  then 

N  Tj 

L(N,  T)  =  CR  +  (N-l)  Cj.  +  C  Z  f  rXt)  dt  (7) 

R  M  r  j=l  0  1 

where  f  =  (Tx,  T2 . TN). 

For  these  purposes,  all  three  corrective  operations  are  assumed 
to  be  instantaneous,  corresponding  to  practical  durations  which  are 
small  compared  to  system  uptime.  The  renewal  cycle  duration  is  then 

N 

D(N,  f)  =  tN  *Z  T.  (8) 

and  the  long-term  average  cost  per  renewal  cycle  is 

C(N,  f)  =  L(N,  f)/D(N,  f)  (9) 
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Theorem  1 


If  i)  r^t)  is  strictly  increasing  in  t,  T1#...,  Tj^ 
ii) 

Then  the  necessary  conditions 


a2r.(t) 

g-rp >0  all  i,  t,  m,  n 
m  n 


(10) 


have  a  solution  corresponding  to  a  global  minimum  of  C(N,  f)  with 
respect  to  the  T^.  Saddle  points,  but  not  a  relative  maximum,  are 
also  possible  solutions  to  (10). 

Proof:  Evaluation  of  the  Hessian  matrix  at  the  stationary  point 

shows  that  a  relative  maximum  cannot  exist  for  any  f. 

C(N,  f )  is  a  continuous  function  of  the  T.  and  it 
approaches  infinity  as  ||f  ||  -*•  0  or  ®.  1 


Theorem  2 

If  the  Tj(t)  have  the  Weibull  form 


i-1 

r,(t)  =  X0[l  +  e  Z 

1  j=l 


Tj]  t® 


(11) 


then  the  optimal  T.  are  ordered  as 

Tf  >  TJ  >  . . .  >  T*  .  (12) 

Proof:  This  result  is  derived  by  showing  that  if  Ti=a<Ti+1=b, 

the  average  cost  will  be  reduced  by  changing  to  Tj  =  b , 

Ti+1  =  a'  w^e  ^eePing  all  other  T,  unchanged.  The 

comparison  is  facilitated  by  the  fact  that  this  inter¬ 
change  does  not  affect  tM  or  r,_(t)  for  k  <  i  or 
k  >  i+1.  N  K 


Another  optimal -solution  property  which  aids  in  computational 


problems  is  summarized  in 


Theorem  3: 


If  a)  r.(t)  are  strictly  increasing  in  t,  and 


b) 


8ri 

at  '  3T 
m  n 


Then  the  optimal  T^  satisfy 
r.(Tp  <  r.(Tp  ; 


all  n  <  i,  m  <  j 


i  <  j 


(13) 


The  Weibull  r.(t)  in  Theorem  2  satisfies  the  conditions  of  Theorem  3. 
The  result  is  exemplified  in  Fig.  1(a)  where  the  peaks  get  higher 
going  to  the  right. 


Proof:  This  follows  from  examination  of  the  following  relation, 

based  on  (10). 

Lt(N,  f)  -  Li+1(N,  f)  =  0  . 

The  preceding  theorems  all  deal  with  the  case  of  fixed  N.  We 
also  want  to  find  the  optimum  number  of  maintenance  operations 
(N-l).  Clearly,  interesting  cases  will  require  CR  >  and  CR  >  Cr- 
The  following  properties  have  been  observed  in  all  numerical  examples 
we  have  studied,  but  sufficient  conditions  and  proofs  have  not  yet 
been  found. 


N 

Conjecture  lm  lN  =  2  Tjf  is  an  increasing  function  of  N . 


Conjecture  2.  The  minimal  cost  as  a  function  of  N  is  convex,  as  in 
Fig.  2,  so  that  an  optimal  N*  can  be  found  by  getting 

f(l)'  f(2) .  until 

C(N+1,  f*N+1))  >  C(N ,  f*N)). 

The  model  considered  here  is  closely  related  to  one  in  Nguyen  and 
Murthy.  They  have  r.(t)  as  functions  of  i  only  (not  on  t.).  In  that 
simpler  case,  equation  (13)  is  replaced  by  r.(Tp  =  r.(Tp  and  Con¬ 
jecture  2  has  been  proved. 

We  now  turn  to  a  slightly  different  problem,  using  (9)  for  the 
criterion,  but  with  case  (b)  of  Fig.  1  to  represent  the  maintenance 
effect.  In  particular 

r^t)  =  X(t  +  t.^)  -  y  ;  0  <  y  <  A  (14) 

Fig.  1(b)  has  y  =  A/2.  Expressing  (14)  in  terms  of  the  intervals  Tj 
we  have 


r.(t)  =  A  t  +  (\-y)  I  T. 
1  1  1 


so  that  r.(t)  is  an  increasing  function  of  t  and  the  T.  (j  <  i). 

Here,  a  necessary  condition  for  the  best  Tj's  based  on  partial 
derivatives  of  (9),  is 

Y(Tj  -  T.)  =  0  ;  i,  j  =  1,2 . N  .  (16) 

Thus,  for  this  model 

Tf  *  Tf  =  ...  T*  =  T*  .  (17) 

Substitution  of  (17)  into  C(N,  T)  leads  to 

T*  *  [-(A-y)(N-1)  +  V(A-Y)^(N-1)"  +  8A(CP  -  (N-DOJ/CJ/2A  (18) 
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In  this  example  we  have  also  found  an  optimal  N  in  every  nu¬ 
merical  example.  Although  the  joint  optimization  of  N  and  f  reduces 
here  to  a  two  parameter  optimization  (N  and  T*),  existence  and 
uniqueness  of  the  overall  minimum  have  not  yet  been  proved. 

All  of  the  preceding  results  can  be  interpreted  with  respect  to 
an  availability  criterion.  If  CR 
times  required,  respectively,  for  replacement,  maintenance  and  minor 
repair,  then  the  long-term  availability  is 

A[N,  f]  = - - =  [1  +  C(N,  f)]"1 

L[N,  t]  +  D[N,  f] 

Thus,  A  is  maximized  here  by  minimizing  C  using  the  techniques 
described  earlier. 

3.  PERIODIC  MAINTENANCE  —  NO  MINOR  REPAIRS 

In  this  model  every  failure  is  followed  by  a  replacement  (re¬ 
newal).  The  replacement  time  has  an  exponential  distribution  with 
mean  value  1/p,  and  maintenance  operations  have  fixed  duration  d. 
The  operating  time  between  maintenances  in  u,  and  failure  may  occur 
during  either  system  operation  or  maintenance.  Failure  rate  reduc¬ 
tion,  of  the  form  of  case  (b)  in  Fig.  1,  occurs  at  the  end  of  each 
maintenance.  The  use  of  equal  operating  intervals  between  mainte¬ 
nances  is  motivated  by  work-rule  convenience  and  the  results  in  the 
previous  section  for  a  related  problem  where  that  structure  is  opti¬ 
mal. 

The  maintenance  effect  is  introduced  differently  here,  but  for 

* 

the  example  of  linearly  increasing  hazard  rate,  there  is  a  direct 


,  CM  and  Cr  are  viewed  as  the  mean 


connection  with  the  models  of  Section  2.  Here  we  say  that  mainte¬ 
nance  reduces  r(t)  by  a  factor: 


r(t p  =  g  r(t[)  ;  0  <  g  <  1  (19) 

Thus,  for  t._1  <  t  <  t.  , 
r(t)  =  rx(t)  -  A. 

^  (20) 
Aj+1  =  (1-g)  [ri(t.)  -  Aj]  ;  =  0  . 

For  the  case  of  r^t)  =  At,  the  maintenance  effect  becomes 


Aj+i  =  (1-g)  [i  As  -  Aj]  (21) 

where  s  =  (u  +  d).  Using  the  notation  of  (4),  this  is  equivalent  to 
r,(t)  =  At  +  8.  ,  ;  i  =  2,3,... 

1  1-1  (22) 
i  i 

0.  =  A  s  I  gJ 
1  j=l 


’  where  the  0j  are  increasing  functions  of  u  (period  between  mainte¬ 
nances). 

We  want  to  find  u  to  maximize  the  probability  that  cycle  avail¬ 
ability  is  above  an  acceptable  level.  With  Tu  and  representing 
the  total  up  and  down  times  in  a  renewal  cycle,  we  have 


A 


c 


u 


VT 


d 


and  the  quantity  tc  be  maximized  is  given  in  (6).  includes  all  of 
the  time  spent  in  maintenance  before  the  failure,  plus  the  replacement 
time  Tr  after  the  failure. 

It  is  convenient  to  define  and  as  the  events  that  failure 

n  n 

occurs,  respectively,  during  the  nth  operating  interval  or  the  nth 
maintenance  interval .  Then 


(23) 


“  Un  Pn(6)  +  Dn  Qn(6) 
P(AC>5]=  I  -  t-  .  *  - ] 


,  U  +  D 

n=l  n  n 


where 


un  *  P'E“1  '  Dn  *  P<En> 

Pn(6)  =  P[AC  >  5|E“}  ,  Qn(5)  =  P[AC  >  6|E*] 


(24) 


Goldstein  (1980)  found  analytical  expressions  for  (23)  in  the  simpler 
case  where  the  durations  of  each  operating  period,  maintenance  period 
and  replacement  were  independent  random  variables  with  exponential 
distributions.  Only  numerical  evaluation  seems  possible  for  the  pre¬ 
sent  model. 

The  piecewise  definition  of  r(t)  in  (20)  can  be  converted  to  a 
corresponding  definition  for  the  lifetime  probability  density  which  can 
be  used  for  numerical  integration  to  evaluate  the  quantities  in  (23) 
for  each  n.  (The  summation  is  truncated  when  terms  become  small.) 

For  t  ,  <  t  <  t  we  have 

n-i  n 


f(t)  =  r(t)  R(t) 
n-1 

R(t)  =  R_(t)  n  R.(t.) 
n  j=l  1  3 


R.(t)  =  exp[-  /  r(x)dx]  :  t.  ,  <  t  <  t. 

J  f  J 

j-1 

Also,  we  note  that  using  TR  for  the  replacement  time  after  failure 

wu 

Pn(5)  =  J  f(t)  P{Tr  <  [t(l-6)/5  -  (n-l)d/6]}dt 
ln-l 


and  a  similar  expression  applies  for  Qn(6). 
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Numerical  examples  have  shown  that  this  approach  is  practical  to 
get  P[AC  >  6]  for  fixed  u,  and  that  a  maximum  can  be  found  by 
varying  the  value  of  u.  Very  small  u  will  introduce  excessive  time 
spent  on  maintenance;  large  u  will  delay  maintenance  until  it  is  too 
late  to  be  effective.  Figure  3  shows  examples  for  P[Ac  >  0.9]  where 
the  optimum  maintenance  period  is  0.07  of  the  mean  lifetime,  or  0.18 
of  the  mean  lifetime,  for  different  parameter  choices. 

4.  CONCLUSIONS 

Optimization  of  maintenance  schedules  has  been  considered  for 
problems  where  maintenance  reduces  the  subsequent  failure  rate  by 
an  amount  depending  on  the  previous  operating  time  of  the  system. 
The  analysis  is  more  difficult  than  in  similar  problems  where  the 
effect  of  maintenance  depends  only  on  the  number  of  previous  mainte¬ 
nances.  Some  useful  properties  of  optimal  solutions  have  been  demon¬ 
strated,  and  others  are  yet  to  be  verified. 
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Case  a 


Numerical  Example,  Convex  C(N,Tj 
Figure  I  Figure  2 

Maintenance-Dependent  Failure  Rotes 


P  [Ac  >  0.9] 


d=O.OI 


d=0.05 


E  [t]  *  8.86  (no maintenance) 

9  *o.rs 


Figure  3-  Optimization  of  Periodic 
Maintenance  Interval  u 


